Analytical description of anomalous diffusion in living cells 
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We propose a stochastic model for intracellular transport processes associated with the activity 
of molecular motors. This out-of-equilibrium model, based on a generalized Langevin equation, 
considers a particle immersed in a viscoelastic environment and simultaneously driven by an exter- 
nal random force that models the motors activity. An analytical expression for the mean square 
displacement is derived, which exhibits a subdiffusive to superdiffusive transition. We show that the 
experimentally accessible statistical properties of the diffusive particle motion can be reproduced by 
this model. 
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The intracellular transport of organelles, vesicles or 
large proteins involves molecular motors that allow the 
fast delivery of cargoes to their correct destination in the 
cell. Molecular motors are proteins able to convert the 
energy from the hydrolysis of ATP in directed motion 
along the cytoskeleton filaments Examples of cy- 
toskeleton motors are kinesin and myosin-V, which move 
along cytoskeleton filaments such as microtubules and F- 
actin 

Single particle tracking techniques have improved sig- 
nificantly in the last years, allowing capturing the 
position of micrometer-sized organelles or beads with 
nanometer and millisecond resolution |3]. Typically, the 
mean square displacement (MSD) of the particle is ana- 
lyzed as a function of the time lag r in order to derive 
the statistical properties of the transport of large car- 
goes within the cell 0, Q , and to analyze the viscoelastic 
properties of the intracellular environment Q . 

Recent experimental works have shown that the MSD 
of the particle exhibit different dynamical regimes on 
different time scales Q- It has been observed that in 
the absence of molecular motors d, 0, [l^, or in the 
case of ATP depletion [ll| the dynamics is subdiffu- 
sive. On the contrary, a crossover from subdiffusion 
(or normal diffusion) to superdiffusion has been reported 
in experiments in which molecular motors are active 

El, [11 [11, Q III, 13 113 ■ In this case, the trans- 
duction of chemical energy into mechanical work pushes 
the cell out of equilibrium 12, [H, [l^, [23| which im- 



plies that the fluctuation-dissipation theorem (FDT) is 
no longer valid. Although it is known that the activity 
of the molecular motors plays a determinant role in the 
observed superdiffusive regime, there is no global model 
accounting for the relationship between the motors ac- 
tivity and superdiffusion up to now. 

In this Letter we propose a stochastic model that takes 
into account the previous facts and enables us to repro- 
duce the main features observed in trajectories of par- 
ticles driven by molecular motors in living cells. For 
this purpose we describe the intracellular transport by 



a generalized Langevin equation (GLE) which includes: 
(i) a delayed friction function that accounts for the vis- 
coelastic properties of medium, (ii) a two terms stochas- 
tic force: a standard internal noise due to thermal ac- 
tivity and an external noise due to active or facilitated 
transport mediated by molecular motors, and (iii) the 
contribution of the experimental errors. We obtain a gen- 
eral expression for the MSD of a particle in a viscoelastic 
environment and in the presence of motor forces which 
can be used to fit experimental data. This approach also 
enables a quantitative description and characterization 
of the different diffusive regimes observed in living cells, 
as was reported in Ref . [l3| . 

The spontaneous motion of a particle immersed in a 
viscoelastic environment is usually described by the gen- 
eralized Langevin equation (GLE) 



mX{t)+ dt' j{t~t')X{t') = F{t) , 



(1) 



where X{t) is the particle position, 7(t) is the dissipative 
memory kernel and F{t) is the random force. 

The integral term accounts for the viscoelastic proper- 
ties of the medium, with the possibility of storing energy 
in the medium and returning it to the particle with a 
finite relaxation time. 

To explicitly include deviation from equilibrium we as- 
sume that the random force F{t) is the sum of two un- 
correlated contributions, i.e. F{t) = ^{t) + x(i), being 
^{t) the standard internal noise due to thermal activity, 
and x{t) an external random force that represents the 
processes that give rise to the active transport. 

The internal noise ^{t), which is responsible for the 
passive motion, is a zero-centered and stationary random 
force with correlation function {^{t)£,{t')) = C{\t — t' \ ) . 
It is related to the memory kernel ^{t) via the FDT [21[ 
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where T is the absolute temperature and ks is the Boltz- 
mann constant. 
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It is now well established that the physical origin of 
anomalous diffusion is related to long-time tail correla- 



tions 
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23l |. In particular, pure power-law correlation 
functions are usually employed to model subdiffusivc pro- 
cess [l^, 24 1 . Then, the noise autocorrelation function 



C{t) can be chosen as 
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where < A < 1, Cq is a proportionality coefficient, tq is 
an arbitrary characteristic time and r(2;) is the Gamma 
function. 

In addition to the thermal noise, we consider an ex- 
ternal contribution originated in the activity of ATP- 
powered motors. This external force x(t), which is not 
related to the dissipation term, is the responsible for the 
FDT violation. In other words, deviation from equilib- 
rium is directly related to the irreversible conversion of 
chemical energy from ATP hydrolysis into t he p article 
motion via the activity of molecular motors [25[. This 
activity, which can be used to generate effectively dif- 
fusive movements by sequences of active directed move- 
ments into random directions, was recently called active 
diffusion (26| and, as we show below, is the origin of the 
transition to a superdiffusive regime. 

Assuming that the network on which the active trans- 
port occurs has a random organization, the random force 
x{t) is chosen as a zero-centered one. On the other 
hand, recent works established that the power spec- 
trum of the noise generated by molecular motors will 



be frequency dependent [11|, [1^, [IJ, |18|. It was also 



established that the autocorrelation function of the to- 
tal noise F{t) has a power- law behavior [ui [H HI . 
Accordingly, we assume a motors force autocorrelation 

A(|i -t'\) = {x(t)x{t')), where A(t) is given by 
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where Aq is a proportionality coefficient and < a < 1. 

It could be thought that the range chosen for a is not 
adequate to reproduce the desired superdiffusive behav- 
ior, and it must be 1 < a < 2. However, in Ref. 27| it 
was established that an external noise with a power-law 
autocorrelation function like (|3|) can lead to a superdif- 
fusive behavior when a is between and 1. This result, 
that has been unnoticed in the literature, will be explic- 
itly shown in this work. Furthermore, considering that 
the power spectrum of the motors force autocorrelation 



is A(w) 



it can be seen that the limit a — > 1 cor- 



responds to a series of instantaneous infinite force pulses 
(white noise limit) while a — > corresponds to the in- 
definitely large memory case, i.e., the so-called strong 
memory limit [28]. Then, an intermediate value of a 
should correspond to a smoothing of discontinuities in 
instantaneous force pulses, as suggested in Refs. 12, T^ . 



This agrees with the well accepted picture of molecular 
motors moving in a step-like manner on microtubules or 
actin filaments [2^ . 

On the other hand, the motion of organelles or vesicles 
is strongly damped in the intracellular media Then, 
the typical damping time constant is too short to be ap- 
preciable experimentally and thus the effect of inertia can 
be neglected in ([T]). In this case, and using the Laplace 
transform technique, the formal expression for the dis- 
placement can be written as 



X{t) ^ xo+ dt'G{t~t')F{t'), 



(5) 



where xq — X{t ~ 0) is the deterministic initial posi- 
tion of the particle. The relaxation function G{t) is the 
Laplace inversion of 

G{s) = , (6) 

where 7(5) = mC{s)/kBT is the Laplace transform of 
the dissipative memory kernel. The relaxation function 
([6|) is independent of the external noise and it is equal 
to the one obtained in the standard internal noise case 
when inertial effects are neglected. 

Typically, the particle trajectory is quantitatively an- 
alyzed in terms of the mean square displacement (MSD), 
which is calculated as {{X{t + r) — X{t))^) where |^(i-f 
t) — X{t)\ is the particle displacement between two time 
points, t denote the absolute time while t is the so-called 
lag time. 

To obtain an analytical expression for the MSD, it is 
necessary to consider the two-time correlation dynamics. 
Using Eq. ([5]) we can write the displacement two-time 
correlation as 

{X{t)X{t')) = xl + I dtiG{t-ti)x 



dt2G{t' -t2){F{ti)F{t2)). (7) 

Since F[t) = ^{t) + x{t), the integral containing the 
correlation function {F{ti)F{t2)) can be split into the 
internal and external contributions. Using relation ([2]), 
and considering the symmetry properties of the correla- 
tion functions C(t) and A(t), the two-time position cor- 
relation function ([7]) can be written as 

{X{t + T)X{t)) = xl + kBT{I{t)+I{t + T)-I{T)) 

+ ( dti{G{ti)H{ti+T) 
Jo 

+ G{ti+T)H{ti)} , 



where 



m 
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Note that, while the relaxation functions G{t) and 
I{t) only depend on the internal thermal noise through 
the memory kernel 7(t), the relaxation function H{t) in- 
cludes the contribution of the external random force. 

For the autocorrelation functions given by ^ and ^ 
the involved relaxation functions can be written as 



where e = Aq/Cq is a dimensionless parameter that mea- 
sures the relative intensity among the motors force and 
the thermal random force. 



m = ^t^tA^ ( - ) , (11) 



Co r(A + 1) \to/ Finally, {(X{t + r) - X{t)f) can be calculated using 

r(f) - ^ /^-LV~^ ® together with Eqs. ^ to Even though the 

To Co r(A) \To/ ' result depends on the relation between A and a, it can 

/ i \ demonstrated that for 2A — a > the MSD have an 

H(t) — eksT —- ( — I , (13) analytical expression given by 

r(A-a-l-l) \tq J 
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= Co/kBT and Bx{a,b) is the incomplete beta 
3l|. While the first term of (|14p represents the 



where 70 
function 

subdiffusive behavior due to thermal activity, the second 
one has its origin on the activity of the external random 
forces. 

Note that the MSD (fT4|) is an aging variable depending 
on the absolute time t and the time lag r i33] . However, 
in typical intracellular tracking experiments an organelle 
or endosome is followed during 10-100 seconds. This time 
is much shorter than the sample preparation durations 
(absolute time). Then, it can be considered that the 
experimental measured MSD is equivalent to the long 
time limit 



MSD{t) = lim {{X{t + r) - X{t)Y) 



(15) 



On the other hand, to make a comparison with exper- 
imental results it is necessary to take into account mea- 
surement errors on the particle position determination 
intrinsic to the SPT experiment or originated in biologi- 
cal activity. It has been established that this effect can be 
introduced by adding an uncorrelated noise of variance 
r]"^ to the mean square displacement [3?| . 

Then, using the asymptotic expansions for the incom- 
plete beta function [SlJ in (jl4p and including the mea- 
surement errors, the MSD (fT5|) can be finally written as 



MSD{t) 
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where 

J., sin(7r(A-a)) - sin(7rA) 
Kx,a = r(a - 2A) , (17) 

TT 

is a positive constant for 2A — a > . 

It is worth pointing out that the second term of (HH) 
is a superdiffusive contribution to the MSD when 1 < 
2A — a < 2. In this case, our model predicts a crossover 
from a subdiffusive to a superdiffusive regimes, with ex- 
ponents A and 2 A — a, respectively. This transition can be 
interpreted as follows: for short enough times the mea- 
surements errors, represented by (2ry)^, dominate, for in- 
termediate time scales a subdiffusive behavior due to the 
viscoelastic properties of the intracellular medium pre- 
vails, while at longer time scales motors activity effects 
dominate leading to a superdiffusive behavior. 

The presented model is characterized by four param- 
eters: A, a, e and 77, where A and a are the power law 
exponents of the internal and external noise correlation 
functions, e is a parameter that measures the relative in- 
tensity between random forces and ry is associated with 
the residual value of the MSD as r ^ 0. Also, the magni- 
tude of the force exerted by the motors can be estimated 
as Fmot ~ \/Ao/r(l — a) where An can be obtained in 
terms of the involved parameters [l3| . 

Interestingly, some recent works have used an empir- 
ical three parameters model of the form A + D*t'^, to 
fit the MSD vs. time lag [sH]- This approach has 
been used indistinctly for systems showing subdiffusive 
{(3 < 1) or superdiffusive (/3 > 1) behaviors. However, 
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as shown above, different regimes can coexist and our 
model allows to describe both situations with a unique 
set of parameters, as we show in Ref.fTl]. For example, if 
in Eq. pB]) the noise (277)^ dominates over the subdiffusive 
term in the measurement temporal range, the empirical 
expression with (3 > 1 holds. On the other hand, the 
subdiffusive behavior observed in the absence of molec- 
ular motors or ATP depletion can be reproduced setting 
e = in (fT6|l . 

In conclusion, we have presented a model that pro- 
vides a physical interpretation of the crossover from 
subdiffusive to superdiffusive behavior observed in single 
particle tracking experiments in living cells. A similar 
approach was recently introduced in the literature 
However, in these works all the forces (internal 
and external) contributions are included in a single term 
and thus, they do not distinguish between the thermal 
and the active forces, a key element to determine 
motor forces in vivo. We believe that the present 
approach can be used to analyze any single particle 
tracking data set obtained in the observation of intracel- 
lular transport driven by molecular motors in living cells. 
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manuscript and helpful discussions. This work was sup- 
ported by grants PICT 928/06 and PICT 31980/05 from 
Agenda Nacional de Promocion Cientffica y Tecnologica, 
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